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ABSTRACT 

We measure the three components of velocity dispersion, aii,ag,a^, for stars within 6 < i? < 30 
kpc of the Milky Way using a new radial velocity sample from the MMT telescope. We combine 
our measurements with previously published data so that we can more finely sample the stellar halo. 
We use a maximum likelihood statistical method for estimating mean velocities, dispersions, and 
covariances assuming only that velocities are normally distributed. The alignment of the velocity 
ellipsoid is consistent with a spherically symmetric gravitational potential. From the spherical Jeans 
equation, the mass of the Milky Way is M (i? < 12 kpc) = 1.3 x 10^^ Mq with an uncertainty of 40%. 
We also find a region of discontinuity, 15 < i? < 25 kpc, where the estimated velocity dispersions and 
anisotropies diverge from their anticipated values, confirming the break observed by others. We argue 
that this break in anisotropy is physically explained by coherent stellar velocity structure in the halo, 
such as the Sgr stream. To significantly improve our understanding of halo kinematics will require 
combining radial velocities with future Gaia proper motions. 

Subject headings: Galaxy: fundamental parameters — Galaxy: halo — Galaxy: structure — Galaxy: 
kinematics and dynamics — Galaxy: formation — stars: Population H 


1. INTRODUCTION 

The global kinematics of halo stars encodes impor¬ 
tant information about the mass, structure, and for¬ 
mation of the Milky Way. The existence of coherent 
struc t ures, such as the Sg r tidal stream (|Ibata et al.l 
11994 iMaiewski et al.ll200^ . demonstrates that at least 
part of the Milky Way halo emerged from the ac¬ 
cretion of smaller galaxie s as expected in hi erarchi- 
cal galaxy formation (e.s.. iSearle fc Zini] 119781 1. The¬ 
oretical simulations suggest that the inner R < 20 
kpc region of the stellar halo should be dominated by 
stars formed in situ whereas the outer region should 


dial orbits (IBullock & Johnstonll2005 

lAbadi et al.l 20061 

Johnston et al. 20081: 

Zolotov et al. 

20091 Font et al.l 

20111: IMcCarthv et al.1120121: IRashkov 

et aLll2013Il. 


Modern studies based on observational data from the 
Sloan Digital Sky Survey (SPSS) support a dual compo¬ 
nent stellar halo (iGaroilo et al.ll20dA 120101: iBeers et al.l 
l2012fl . In this view, the inner halo i? < 15 kpc exhibits 
a flattened distribution of stars on radial orbits with no 
net rotation contrasted with the outer halo i? > 20 kpc 
characterized by a sp herical distribution of stars w ith net 
retrograde rotation. iSchonrich et al.l (|201lL [20l3l . how¬ 
ever, argue that the dual halo results from observational 
errors and selection effects, and cannot be distinguished 
from a single halo full of substructure. 

Intriguingly, many other observers find a discontinuity 
in both the number density and the orbital properties 
of halo stars near R 20 kpc. Star counts of main se - 


quence turn-off s tars (iBell et al 


RR Lyrae stars (jWatkins et al.ll2009y Sesaj_e1^_aj^ 
blue horizont al branch (BHB) s tars (|Deason et ah 


Sesar et al.l 


2011 ) 


2013) 


Will 


and K giants (iKafle et al.ll20l4 l all exhibit a break in the 
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number density around R = 16 — 26 kpc. Radial velocity 
surveys of BHB stars impl y systematically rad i al orb its in 
the inner and outer halo (|Deason et al.l[2MTL l2012h but 
more tangen^tial orbits in the region 15 < i? < 25 kpc 
(jKafle et al.ll2012l : iDeason et al.ll201^ . The anisotropy, 
13, that depends on the ratio of tangential and radial 
velocity dispersions, provides a useful means for quanti¬ 
fying systematic velocity changes in the break region. 

Several studies have exploited large samples of halo 
stars to measure the velocity dispersion and anisotropy 
profiles be y ond R > 10 kpc with varying results. 
iSirko et ^ (|2nn4f ) employed a sample of 1,170 BHB 
stars selected from the SDSS Data Release 4 to mea¬ 
sure the anisotropy [3 = 0.1±0.2, a result consistent with 
isotropy, for stars with a median distance from the Galac¬ 
tic center of i? ^ 25 kpc. IDeason et al.1 (|2012t) analyzed 
1,933 BHB stars from SDSS Data Release 8 (DR8) with 
16 < i? < 48 kpc to And a radially biased anisotropy of 
(3 = 0.51q'2®. Looking for evidence of a multi-component 
halo, IKafle et ^ (I2012D analyzed an SDSS DR8 sam¬ 
ple of about 4, 500 BHB stars and found the anisotropy 
was radially biased, (3 = 0.5, for 9 < i? < 12 kpc and 
25 < i? < 56 kpc. To their astonishment, they discov¬ 
ered a sharp dip, j3 ~ —1.2, in the anisotropy parameter 
profile at i? ~ 17 kpc that they could not explain as 
arising either from halo substructures or from accretion. 

Stellar velocities also provide information about the 
gravitati onal potentia l and mass of the Galaxy. The 
spherical iJeansI (jl915fl equation provides a quantitative 
link between observ ations and the underlying gravita¬ 
tional potential ('e.g.. lBinnev fc Tremainell2008D . Several 
sophisticated Jean analyses have been performed in re¬ 
cent years, exploring the importance of density profiles, 
anisotropy as sumptions, and potential models with SPSS 
observations llDeason et ani2012l : IKafle et ai]l2012ll20l4 

iLoebman et al.ll2014fl . This paper explores the kinemat¬ 

ics of stars with 6 < i? < 30 kpc using new and existing 
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and SDSS (gray) surveys in Galactic latitude and longitude. The 
SDSS BHB stars have essentially the same footprint as the SDSS 
F stars (gray). 

data sets. 

We present a spectroscopic radial velocity survey of 
6,174 faint 18 < r < 21 F-type stars obtained with 
the Hectospec spectrograph on the 6.5m MMT telescope. 
We target F-type stars because they are the densest lu¬ 
minous stellar tracer at heliocentric distances of 12 to 
20 kpc. Padova tracks predict that a 10 Gyr old 0.8 
Mq star with [ Fe/H] = —1.7 has an absolute magnitude 
of Mr = +4.5 (IGirardi et al.l I2002t IMarigo et al’l 120081 
iBressan et al.l 120121) . a value consisten t with globular 
cluster observations (|Newbv et al.ll201ll) . Thus r = 20 
and 21 mag F-type stars probe the halo at 12 and 20 
kpc heliocentric distances, respectively. To validate our 
results and improve our statistics, we also make use of 
13,480 F-type stars culled from the SDSS. Because the 
SDSS F star sample is shallower than our Hectospec F 
star sample, we also use 3,330 BHB stars culled from the 
SDSS to better constrain the more distant region i? > 15 
kpc. 

In the next section, we describe the three data sets 
used in our analyses. Section 3 discusses the geometry of 
projecting the components of velocity on the line of sight 
and develops maximum likelihood statistical methods to 
infer the distribution of velocities from the observed line- 
of-sight velocities. Section 4 describes the mean velocity, 
dispersion, covariance, and anisotropy profiles; the tilt 
of the velocity ellipsoid; and the mass of the Galaxy. In 
Section 5, we test the assumption of independent, nor¬ 
mally distributed velocities and investigate the potential 
effect of the non-normal velocity distribution of the Sgr 
stream on our analysis. We summarize our conclusions 
in Section 6. 


2. DATA 

We analyze three data sets, including a new sample 
of radial velocities obtained with the Hectospec spec¬ 
trograph on the 6.5m MMT telescope. To validate our 
results and improve our statistics, we add F-type stars 
culled from the SPSS Stellar Parameter P ipeline (SSPP, 
iLee et al']l2008 : lAllende Prieto et al.ll2008[) and the BHB 
star sample of IXue et al.l ( 201111 . 


2.1. Hectospec F Star Sample 
2.1.1. Observations 

Hectospec is a 300 fiber, multi-object spectrograph 
with a 1 d egree diameter field of view on the 6.5m MMT 
telescope (iFabricant et al.l [20051) . All observations are 



X (kpc) 


Fig. 2.— Spatial distribution of F stars in the Hectospec (blue) 
and SDSS (gray) surveys projected onto the Galactic Cartesian 
coordinate X-Z and X-Y planes. The Sun is located at X = —8 
kpc. 


made with the 270 line mm“^ grating, which provides a 
spectral resolution of 5 A over the spectral range 3700 - 
9100 A. Hectospec fibers are assigned such that high pri¬ 
ority targets are assigned first, followed by lower priority 
targets. All targets and their prio rities are determined 
on the basis of SDSS de-reddened (iSchlegel et al.l [19^ 
point spread function magnitudes and colors. We indi¬ 
cate de-reddened magnitudes and colors with a subscript 
0 . 

We acquired 3,197 spectra between 2004 April and 
2005 July, the first year of Hectospec operations, with 
a dedicated halo star observing program. Our dedicated 
program targeted 17 < 50 < 20 mag stars, prioritizing 
A-type (candidate BHB) stars and filling the remain¬ 
ing Hectospec fibers with increasingly red stars out to 
(5 ~ ''’)o < 0.5. Given the relative surface densities of 
stars, 80% of the targets were F-type stars. 

We acquired an additional 8,143 spectra between 2009 
January and 2015 July with a parallel Hectospec observ¬ 
ing program. Our parallel program took advantage of all 
unassigned fibers in Smithsonian Astrophysical Observa¬ 
tory (SAO) Hectospec programs, and filled those fibers 
with 18 < ro < 21 F-type stars. Our target selection 
for the parallel program prioritized the faintest stars at 
the main sequence turn-off {g — r)o — 0.25. The parallel 
program yielded between 200 and 2,000 spectra a year, 
depending on fiber availability and Hectospec usage. 
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Fig. 3. — Color-magnitude distribution of the Hectospec sam¬ 
ple of 6,174 F-type stars. All photometry comes from SDSS de- 
reddened point spread function magnitudes. The median depth is 
ro = 19.0, or 8 kpc. Applying a |Z| > 5 kpc cut removes the stars 
approximately below the dashed line. 

Given the nature of our parallel observations, the over¬ 
all angular distribution of F stars is irregular and clumpy, 
as seen in Figures [T] and [H The inner stellar halo of the 
Milk y Wav is expected to be sp atially well-mixed, how¬ 
ever (jBullock fc .Tohnstonl 120051) . Thus the F stars, se¬ 
lected by color, should be fair probes of halo kinematics. 

2.1.2. Radial Velocities 

All Hectospec spectra are processed by th e SAP Tele¬ 
scope Data Center’s data reduction pipeline (jMink et al.l 
120071) . We visually inspect each spectrum for quality 
control. The major contaminants are M-type stars (4% 
of the spectra) and quasars and miscellaneous galaxies 
(1% of the spectra). All contaminants are removed from 
the final Hectospec sample. 

We measure stellar radi al velocities with the RVSAO 
cross correlation package (jKurtz fc MinklflO^ . We use 
cross-correlation templates constructed from Hectospec 
observations of 33 bright A- and F-type stars with known 
velocities from the Century Surv ey Galactic Halo Project 
(jBrown et al.ll2003Ll2005bl l2008D . We adopt the weighted 
mean velocity for objects with more than one observation 
(7% of the spectra). Our median radial velocity uncer¬ 
tainty is ±15 km s“^ for tq = 20 mag F-type stars. 

A cross-identification search finds 566 spectra in com¬ 
mon with the SDSS SSPP catalog. The mean differ¬ 
ence between our radial velocity and SSPP elodiervfinal 
is 1.0 ± 20.4 km s“^. Our systematic error with respect 
to SDSS is thus 1 km s“^, which is consistent with the 
sum of the measurement uncertainties. 

2.1.3. Final Sample 

Our targets were selected from a variety of SDSS data 
releases. We unify the photometry for this paper using 
SDSS Data Release 10 de-reddened point spread func¬ 
tion magnitudes and colors. Figure |3] presents a color- 
magnitude plot of the cleaned Hectospec sample, which 
has a median depth of tq = 19.0 mag. 


Our spectra were acquired with a variety of exposure 
times and observing conditions. We therefore remove low 
quality observations as the first step of making a clean 
sample. We require that each spectrum has more than 60 
instrumental counts and a signal-to-noise ratio per pixel 
(S/N) > 4 in the continuum at 5000 A. Radial velocity 
measurements are possible from lower quality spectra, 
but our S/N requirement ensures repeatable velocities 
with <40 km s“^ errors. Second, we require that the 
ratio of counts in the continuum at 7800 A and 5000 A 
is < 1.25. This requirement, given the Hectospec grating 
blaze function, eliminates all non-stellar objects, such as 
quasars, and spurious red objects, such as M-type stars. 
Finally, we impose 0.15 < {g — r)o < 0.45 and 0.4 < 
{u — g)o < 1.4 to yield a clean sample of 6,174 F-type 
stars. 

Given the low spectral resolution and modest S/N ratio 
of the Hectospec spectra, we do not p erform stellar atmo - 
sphere model fits. Instead, we use the llvezic et al.l (l2008f ) 
photometric paral lax relation in combination with the 
IBond et al.l (|2010ll metallicity calibration. These rela¬ 
tions use de-reddened broadband ugri colors to estimate 
stellar luminosity Mr and metallicity [Fe/H]. Our radial 
velocity cross-correlation templates span A through F 
spectral types, and the best-matching template to each 
spectrum validates this approach. 

The average [Fe/H] = —1.4 and Mr = ±4.5 for our 
entire sample are consistent wi t h expectations of halo 
F-type stars (llvezic et al.l 1200^ IBond et ^ 120101) . At 
faint r > 19 magnitudes and low S/N ratio, however, 
the scatter in our estimated [Fe/H] becomes implausibly 
large. Spectroscopic metallicity measurements show that 
F-type stars at these dept hs have —3 < [Fe/H] < — 1 
and mean [Fe/H] = —1.6 (jAllende Prieto et al. 120141) . 
To minimize the distance error for stars with poorly es¬ 
timated [Fe/H], we assign [Fe/H] = —1.6 to stars with 
[Fe/H] < — 3 or [Fe/H] > — 1. We expect t hat our dis¬ 
tance estimates me pr ecise to about 15% (jlvezic et al.l 
120081: IBond et "ni l2010fl . Our measured heliocentric ra¬ 
dial velocities and estimated distances are provided in 
the Appendix lA.31 

To remove all significant disk contamination, we im¬ 
pose \Z\ > 5 kpc. This cut is motivated by the observed 
distributions of velocity and metallicity for F stars, which 
show that the dis k population becomes negligible around 
\Z\ =4 — 5 kpc (llvezic et al.l 120081: iGarollo et al.ll2?)Tol : 
IBond et all 120101 : iFermani fc Schonricbll2013D . Restrict- 
ing the sample to stars with Galactic rest frame velocities 
\vrf\ < 500 km s ~^ to eliminate pote ntial unbound hy¬ 
pervelocity stars (jBrown et al.l l2005al) does not remove 
any stars from the sample. Our final Hectospec sample 
contains 3,049 F-type stars. 

2 .2. SDSS F Star Sample 

The SSPP catalog provides spectroscopic measures 
of effective temperature, surface gra vity, and metal- 
licity for stars observed by SDSS ([Lee et al.l 120081 : 
lAllende Prieto et'aI1l2008[) . We construct an F star sam¬ 
ple from the SSPP catalog by selecting stars that have a 
stellar classification of F and extinction corrected colors 
satisfying 0.4 < {u — g)o < 1.4, 0.15 < {g — r)o < 0.45, 
and —0.2 < (r — i)o < 0.6. We esti mate heliocentric dis¬ 
tances using the llvezic et al.l (|2008D photometric parallax 
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relation for consistency with the Hectospec sample. The 
median depth of the resulting F star sample is rp = 17.4, 
or d ~ 4 kpc. 

Requiring \vrf\ < 500 km s“^ to avoid potential con¬ 
tamination by hypervelocity stars removes five stars, all 
of which have large radial velocity errors. Removing all 
objects in common with the clean Hectospec F-star sam¬ 
ple and performing a cut \Z\> h kpc to remove potential 
disk contaminants yields a sample of 13,480 stars. The 
median heliocentric radial velocity uncertainty is ±11 
km s“^ for r = 20 mag F stars. 

Figured] compares the angular distribution of the Hec¬ 
tospec and SDSS F star samples on the sky. The overall 
distribution reflects the SDSS imaging footprint, which 
covers a large range of Galactic longitude but a restricted 
range of Galactic latitude, convolved with the spectro¬ 
scopic survey regions. Figure [2] displays the correspond¬ 
ing spatial distribution of these stars projected onto X-Y 
and X-Z planes in Galactic Gartesian coordinates. The 
majority of F stars are within about 10 kpc of the Sun 
and densely probe out to i? ~ 20 kpc Galactocentric 
distances. 

2.3. SDSS BHB Star Sample 

IXue et aP (|2011ll spectroscopically identify 4,985 stars 
in SDSS as luminous BHB stars. These objects are 
evolved, metal-poor halo stars with typical absolute mag¬ 
nitudes of Mg ~ ±0.8 mag, significantly more luminous 
than the F-type stars. Because the BHB stars come from 
SDSS spectroscopic fields, the sky coverage is essentially 
identical to that of the SDSS F star sample (Figure dj)- 

We esti mate BHB ab s olute magnitudes and distances 
using the iDeason et al.l (j2011f ) color-magnitude relation. 
The median apparent magnitude of the BHB sample is 
go = 16.7 mag, corresponding to a depth of 15 kpc. Im¬ 
posing a restriction on velocity that \vrf\ < 500 km s“^ 
to avoid potential contamination by hypervelocity stars 
does not eliminate any stars. Retaining only stars with 
|Z| > 5 kpc to avoid disk contamination yields a sample 
of 3,330 stars. The median heliocentric radial velocity 
uncertainty is ±9 km s“^ for g = 19 mag BHB stars. 

2.4. Galactic Rest Frame Velocity 

We transform all heliocentric velocities, Vheiio, to 
Galactocentric rest frame velocities, Vrf, assuming a 
circular velocity o f 235 km s~^ (IReid et al.l 120091: 
McMillan fc Binned 120101 : iBovv et al.ll2012t IReid et^ 
2014H an d a solar motion of (U V, W) = (11.1,12.24, 7.5) 
km s“^ (|Sch6nrich et al.l[2010l) . 

Vrf = Vfieiio ± 235 sin(l) cos(6) ± 11.1 cos(l) cos(6) 

± 12.24sin(l) cos(6) ± 7.25sin(6). (1) 

Figured] plots summary statistics for the observed Vrf 
distribution for the combined data sets. We use bins of 
fixed 4 kpc width and plot points at the mean R of each 
bin. Figured] (panel a) shows that the velocity dispersion 
of the combined data sets is relatively constant except for 
the range 15 < i? < 18 kpc, where the velocity disper¬ 
sion declines suddenly from 123 km s“^ to 102 km s“^. 
Interestingly, the Hypervelocity Star survey measures a 
113.9 ± 6.6 km s~^ velocity dispersion 15 < i? < 20 kpc 
(|Brown et al.l 1201^ in perfect agreement with the aver¬ 
age velocity dispersion in this region. 




mean R [kpc] mean R [kpc] 




Fig. 4. — Statistics of the observed v^f distribution for the com¬ 
bined sample versus the mean R for each bin. Bins have fixed 4 
kpc width and thus overlap in R, with varying number of stars 
per bin. Panel a) is the velocity dispersion, panel b) is the mean 
velocity, panel c) is the kurtosis, and panel d) is the p value of 
Kolmogorov-Smirnoff (K-S) and Anderson-Darling (A-D) tests for 
normality of the v^-f distributions in each bin. 

Figure d] (panel b) shows some apparent variation in 
the mean velocity of the stars, but for most bins the 
mean velocity does not differ significantly from 0 km s“^. 
Figured] (panel c) shows that the kurtosis of the distri¬ 
bution differs significantly from zero, however, indicating 
that the distribution of Vrf departs from a normal distri¬ 
bution for bins i? < 18 kpc. Both Kolmogorov-Smirnoff 
and Anderson-Darling tests significantly reject a normal 
distribution for bins i? < 18 kpc. We caution that both 
statistical tests are sensitive to the number of stars and 
note that the outer bins contain many fewer stars than 
the inner bins. 

3. THEORY AND METHODOLOGY 

There are few tests of the dynamical models used to 
derive Galactic parameters. Halo models typically pos¬ 
tulate a high degree of spatial symmetry and assume that 
the velocities of halo stars are normally distributed and 
uncorrelated with zero means. The presence of structure 
and star streams expected in hierarchical galaxy forma¬ 
tion might violate the standard assumptions. To explore 
these ideas, we calculate the velocity means, dispersions, 
covariances, and the anisotropy of the inner halo from 
observational data using a minimal set of assumptions. 

Ideally, our input would be 3-dimensional velocities 
for a large number of stars covering large areas of sky. 
Unfortunately, t he typical ±5 ma s yr“^ proper mo¬ 
tion uncertainty (|Monet et al.ll200§] of an r = 20 mag, 
d = 12 kpc F-type star translates into a ±284 km s“^ 
uncertainty in tangential velocity. Tangential velocities 
are consequently uninformative except in large statisti- 
cal averages, which remain sensitive to systematic error 
(|Fermani fc Schonrichl[20T^ . Radial velocities are about 
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Fig. 5. — Viewing geometry depicting the relation between he¬ 
liocentric and Galactocentric coordinates systems. Gray spherical 
caps show i? = 20 kpc for \Z\ > 5 kpc. Solid circles show where 
stars with heliocentric distances of d = 15 (red), 20 (purple), and 
25 (blue) kpc fall on this i? = 20 kpc surface. 

25 times more accur ate. The typical ±11 km s“^ radial 
velocity uncertainty (|Lee et al.ll2008fl for an r = 20 mag 
star makes radial velocity our tool of choice. Because 
the Sun lies about 8 kpc from the Galactic center, we 
use the heliocentric radial velocity, combined with the 
angular position and the distance of each star, to con¬ 
strain Galactocentric tangential velocity components. 

3.1. Projection Factors 

Geometrically, the observed line-of-sight velocity, vios, 
is the projection of its Galactic velocity components, 
V ={vr,vo,v^), on the line-of-sight unit vector, n, 

Uios=v-n= ^ ViPi, (2) 

where pi are projection factors for each component. We 
use Galactic spherical coordinates because we are dealing 
with the halo. The projection factors depend on angular 
position {l,b), heliocentric distance d, and the distance 
from the Galactic center R (see Appendix lA.ll for de¬ 
tails). Figure [5] illustrates the geometry. Sight lines from 
the Sun at different Galactic longitude, Z, and latitude, 
6 , intersect the gray R — 20 kpc surface at points cor¬ 
responding to different distances from the Sun, although 
they are all the same distance from the Galactic center. 

Importantly, different directions on the sky provide in¬ 
formation about different Galactic velocity components. 
This point is widely understood but rarely quantified. 
Gonsider the region on the sky |6| > 80°: 70% of vios 
is in the 9 component for stars at heliocentric distance 
d = 5 kpc, but only 20% of vios is in the 6 component for 
stars at d = 15 kpc; the remainder is in the R component. 
Recovering tangential velocities from vios clearly requires 
large observational samples with broad sky coverage. 

Figure [6] quantifies the projection factors from the ob¬ 
server’s perspective as functions of angular position (Z, b) 
and depth. The three columns in FigurelHlare for the pn, 
pe, or p^ projection factors. Golor indicates the mag¬ 
nitude of the projection factor. The white contours in 


each panel mark the regions on the sky at heliocentric 
distances of d = 15, 20, and 25 kpc. The three rows in 
Figure [6] represent Galactocentric distances R = 10,15, 
and 20 kpc. The sky coverage reflects the jZj > 5 re¬ 
striction as R increases. 

As expected, information about vg comes primarily 
from the Galactic polar regions; information about v^j, 
comes primarily from regions near Galactic latitudes of 
60° and 300°. At large distances, the observed line-of- 
sight velocity contains little information about the tan¬ 
gential velocity components and becomes essentially ra¬ 
dial. If there is any correlation between angular position 
and kinematics, a systematic bias is introduced. 

3.2. Maximum Likelihood Estimation 

We estimate the velocity means, dispersions, covari¬ 
ances, and anisotropy of stars in the halo as functions of 
distance from the Galactic center using a maximum like¬ 
lihood procedure. We adopt the standard assumption 
that the halo star velocit ies have norm al di stributions 
(|Binnev fc Tremainell2008ll : see ADDendix lA.2l for details. 
This minimal assumption allows us to test whether or not 
the velocity components are normally distributed, have 
zero means, and are uncorrelated. Correlations among 
velocity components, if they exist, may reveal underly¬ 
ing dynamic structures in the halo, such as star streams. 

Standard practice is to estimate the velocity disper¬ 
sions by binning observations in a series of contiguous 
intervals in R. This procedure yields a series of discrete 
estimates, one for each bin. We determine model param¬ 
eters by maximizing the log likelihood function using an¬ 
alytic firs t and second derivatives and the R sta tistical 
software (|R Core Teaml[20H : iToomet et al.ll2012ll . 

We validate our maximum likelihood calculations 
against a simulated set of stars with 6 < i? < 30 kpc and 
normal velocity distributions. Our maximum likelihood 
calculation finds the correct solution for the underlying 
model parameters. 

We estimate both a restricted and a full statistical 
model. In the restricted model, the velocities are as¬ 
sumed to have zero means and to be uncorrelated, i.e., 
the covariance matrix is diagonal, as in previous stud¬ 
ies of velocity dispersion in the halo. Estimating the 
full statistical model allows us to test empirically the as¬ 
sumptions of the restricted model that the velocities have 
zero means and are uncorrelated and thus test the funda¬ 
mental assumptions underlying dynamical models of the 
Milky Way halo. For example, if the gravitational po¬ 
tential is spherical ly symmetric, the velocity -dispersion 
tensor is diagonal (iBinnev fc Tremain3l2008l l. 

4. RESULTS 

We begin by investigating the velocity dispersion and 
anisotropy of the Milky Way halo using our independent 
Hectospec F star sample. We then combine our F star 
sample with the SDSS data sets to test our results with 
a larger set of observations. 

4.1. Anisotropy Parameter 

One way to characterize the orbital structure of a 
spherical system, such as the h alo, is through the 
anisotropy parameter defined as (IBinnev fc Tremaind 
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Fig. 6.— Line-of-sight projection factors at different Galactocentric distances in the halo. Columns show the components pg, and 
P(k. Rows show the projection factor values at R = 10, 15, and 20 kpc. For reference, white contours mark heliocentric distances of d = 15, 
20, and 25 kpc. 


[2tM eq. 4.61) 


/3 = 1- 


2^1 


(3) 


Because the anisotropy parameter depends on the ratio of 
the velocity dispersions squared, it rapidly becomes more 
negative as the tangential velocity dispersions increase 
relative to an- The value of the anisotropy parameter is 
/3 = 1 for perfectly radial orbits, /3 = 0 for isotropy, and 
P = —oo for perfectly circular orbits. From our velocity 
dispersions estimates, we can quantify the anisotropy of 
the orbits in the Milky Way halo and investigate the 
break in 0 found in previous at stud ies at i? ~ 20 kpc 
(jKafle et al.ll2012l : iDeason et al.ll201^ . 


4.2. Results from Hectospec F Star Sample 

Using the line-of-sight velocities and positions from our 
Hectospec F star sample, we explore the Galactic ra¬ 
dial profile of the velocity dispersions cr/j, erg, cr$, and the 
anisotropy parameter, 0. We begin with the restricted 
model, which assumes that the velocity components have 
zero means, Hi = 0, and are uncorrelated, = 0 for 
i ^ j- We relax both of these assumptions when we 
investigate the larger combined data set in Section 

To estimate the parameters of our statistical model, 
we partition our observations into four contiguous radial 
bins containing ribin = 762 or 763 stars and estimate the 
parameters and their errors separately for each bin us¬ 
ing maximum likelihood. A caret symbol distinguishes 
our parameter estimates from the parameters themselves. 
Figure 0 and Appendix IA.4I present our estimates of the 
velocity dispersions and the anisotropy parameter. The 
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Fig. 7.— Velocity dispersions and anisotropy as function of R for 
the Hectospec sample. The 3,049 stars in the Hectospec sample are 
partitioned into four radial bins of equal numbers. The estimates 
are plotted at the mean value of R for stars in each bin. Horizontal 
error bars show bin size; vertical error bars show Icr uncertainties. 
Appendix I A. 41 presents the numbers. 
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Fig. 8.— Mean velocity components of the combined sample. 
The 19, 859 stars in the combined sample are partitioned into eight 
radial bins with equal numbers of stars (green circles) and also into 
overlapping bins of fixed 4 kpc width (red triangles). Horizontal 
error bars show bin size; vertical error bars show la uncertainties. 
The estimate of jiff, drops significantly between 15 < < 25 kpc. 


Fig. 10.— Velocity dispersions and anisotropy for the combined 
sample. Sample, error bars, and bins as in Figure [8] The fixed 
width bins (red triangles) overlap in R and have a varying num¬ 
ber of stars per bin, whereas the fixed number bins (green circles) 
have varying widths that do not overlap. Either way, there is a 
large increase in estimated ag and (T^, and a corresponding drop 
in anisotropy /3, in the range 15 < < 25 kpc. 
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Fig. 9.— Velocity covariances for the combined sample. Sample, 
error bars, and bins as in Figure [8] Note the significant positive 
covariance between vji and vg in the region /? = 14 — 17 kpc. 

solid points in Figure [7] represent the value of the point 
estimate for the observations in a bin plotted at the mean 
R for stars in that bin. The vertical lines mark the 
Icr confidence intervals for the estimates. Note that the 
width of the bins varies as a function of Galactic radius, 
i?, due to the declining number density of stars and the 
limiting magnitude of the observations. Taken together 
the graphs reveal how the velocity dispersions and the 


Fig. 11.— Contours of constant illustrating the correlation 
between the velocity dispersion components in the 14 < it < 18 
kpc bin. The dispersion shows the largest range of values and 
thus is the most poorly constrained, while an and ag are better 
constrained but correlated. 

anisotropy parameter vary with Galactic radius for the 
F stars in our sample lying more than 5 kpc above the 
Galactic plane with 6 < i? < 30 kpc. 

The results for an and those for the first two points 
of and /3 are consistent with previous findings. 

The drop in the anisotropy parameter occurs near the 
discontinuity observed previously by others. Large values 
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for the tangential velocity dispersions ag and lead to a 
negative value for /3 implying that orbits are tangentially 
biased at these radii. 

These apparent anomalies could arise from several 
sources. The stellar velocities may not be normally dis¬ 
tributed, or they may be correlated. The presence of 
structure or star streams in the inner halo may account 
for the discrepancies by violating the underlying assump¬ 
tions of the spatial symmetry and statistical distribution 
of stellar velocities. 

4.3. Results from Combined Sample 

Combining the observations of the Hectospec and 
SDSS samples provides greater statistical power. With 
the larger number of observations, we can explore a more 
complete model of the statistical distribution of halo 
stars. We relax the assumptions of the restricted model 
and include estimates of the means of the spatial veloci¬ 
ties, p,i, to address the possibility of systematic motions, 
and the off-diagonal elements of the variance-covariance 
matrix, to allow for correlated velocities. 

Figures [8] - [10] and Appendix IA.4I present the max¬ 
imum likelihood estimates of the mean velocities, ve¬ 
locity dispersions, covariances among the velocities, and 
anisotropy parameter of the stars in the combined sam¬ 
ple. The combined sample includes 19, 859 stars span¬ 
ning 6 < i? < 30 kpc. We estimate the parameters as a 
function of Galactic radius using two methods for parti¬ 
tioning the data: 1) placing approximately equal num¬ 
bers of stars (ribin = 2,482 or 2,483) within each interval 
in R (green circles) and 2) using bins with a fixed width 
of 4 kpc (red triangles). The error bars represent Icr. 
The two partition methods yield similar results. 

For i? < 15 kpc, the estimated mean velocities are 
small and, given the estimated errors and the uncertain¬ 
ties in the solar motion, are consistent with zero (Figure 
0 ). Within the range 15 < i? < 25 kpc, however, and 
fig drop significantly below zero. Within 15 < i? < 25 
kpc, the stars have mean velocities fi^ ~ —50 to —80 
km s“^ and fig ~ —15 to —40 km s“^. 

The estimated off-diagonal covariances among the ve¬ 
locity components (Figure |9|) are consistent with zero 
at the 2(7 level with some exceptions. Positive covari¬ 
ances Tirq between vr and vg and between vr and 
occur in the region i? = 14—17 kpc, for example. 
The covariance estimates typically are at least an order 
of magnitude smaller than the variances, except for the 
covariance Yig^ between vg and v^f, beyond i? > 20 kpc 
(Figure |9|), which is poorly estimated. The positive cor¬ 
relation between vg and occurs in approximately the 
same region where the mean velocity, /i^, turns signifi¬ 
cantly negative, again suggesting the possibility of cor¬ 
related motions in this region. 

The velocity dispersion estimates exhibit unexpected 
behavior in the range 15 ;< i? < 25 kpc. Over the entire 
span of 6 < i? < 30 kpc, the radial velocity dispersion, 
aR, falls from about 150 km s“^ to a low of 50 km s“^ and 
then recovers to about 90 km s“^. The changes in the 
estimates dg and are more dramatic. Over the entire 
span, ag jumps from about 110 km s“^ to a maximum of 
about 360 km s“^; increases from about 110 km s“^ 
to a peak of about 250 km s“^. Similarly, the anisotropy 
estimate, /3, declines from 0.5 to around —20. 



Fig. 12.— Comparison of anisotropy results. We plot the loga¬ 
rithm of (1 — /3) to present a more balanced comparison of tangen¬ 
tial and radial anisotropies. We present anisotropy estimates for 
our overlapping fixed 4 kpc bins (red squares), together with the 
anisotropy estimates of prior researchers. 

To illustrate the correlations among the velocity dis¬ 
persion components. Figure [TT] plots contours of constant 
for stars in the 14 < i? < 18 kpc bin. The velocity 
dispersion estimate (t^ shows the largest range of values 
and thus is the most poorly constrained of the three ve¬ 
locity components. The dispersions or and ag are better 
constrained but correlated. Numerically, the maximum 
likelihood calculation allows the two tangential compo¬ 
nents to settle around 300 km s“^ as (tr goes to zero. 
This is not a physically plausible scenario. The high 
velocity dispersions for ag and a^ and the positive co- 
variance in the range 14 < i? < 18 kpc imply that some 
of the stars must have velocities exceeding the Galactic 
escape velocity. 

We compare our results for the anisotropy with previ¬ 
ous research in Figure [T2| We plot our /3 estimated in 
fixed width bins of 4 kpc. Our results at i? < 15 kpc and 
at 77 > 25 kpc are i n good agreement with prior research 
(iSirko et al.l 2004 iBond et al.l 120101 : iKafle et all 120121 : 
iDeason et al. I20l'^ 1201311 ^ In the range 15 < i? < 25 
kpc, however, our estimates of /3 are substantially more 
negative (corresponding to tangentially biased orbits). 

We attribute the less negative values of /3 estimated 
by other researchers to their use of much wider bins 
in R, which smooth the actual dispersion profiles; to a 
lower Z cutoff of 4 kpc, which may increase contami¬ 
nation by disk stars; and to different statistical mod¬ 
els that marginalize over tangential velocities and, in 
some cases, also impose additional assumptions about 
the gravitational potential and number density of stars. 
As shown by the hori z ontal bars in R , the e s timates of 0 
bv ISirko et al.l (1200411 , IDeason et all (I2012D . iKafle et al] 
(|2012f) . and IDeason et al.l (|2013f ) are based on broad 
ranges in R that extend beyond the interval where our 

























































Stellar Velocity Dispersion of Inner Halo 


9 


results diverge. Although iKafle et al.l (|2012ll do not re¬ 
port the sizes of the bins used for their estimations, we 
deduce from their original data that their bins increase 
in width from approximately 5 kpc at i? = 17 kpc to 7 
kpc at i? = 23 kpc. 

At first glance, our results corroborate the discontinu¬ 
ity in /3 observed around i? ~ 20 kpc. Our results suggest 
the presence of correlated stellar motions within the re¬ 
gion of 15 < i? < 25 kpc, perhaps resulting from star 
streams or other structure. The implausibly high veloc¬ 
ity dispersions also suggest that there may be another 
explanation for the discrepancy, namely that our under¬ 
lying statistical model and those of previous researchers 
may be misspecified. We investigate these possibilities 
further in Section 5 after first considering the robustness 
of our results, assessing the alignment of the velocity el¬ 
lipsoid in the halo, and estimating the mass of the Milky 
Way interior to i? = 12 kpc. 


4.4. Robustness of Results 

Our results are subject to various sources of poten¬ 
tial error. We test the robustness of our parameter esti¬ 
mates with a series of sensitivity analyses on the input 
data. Systematically overestimating distances, for ex¬ 
ample, may inflate estimates of the tangential velocities 
at larger distances causing larger estimates of the tan¬ 
gential velocity dispersions, ve and v^, and negatively 
biasing estimates of the anisotropy (3. To investigate the 
effect of systematic errors in stellar distance estimates, 
we recompute our results with all distances increased and 
decreased by 20%. To investigate the effects of our choice 
of the location of the Sun with respect to the Galactic 
center, we recompute our results for Xq = —8.5 kpc. To 
investigate our choice of the Sun’s circular velocity, we 
recompute our results for circular velocities of 200 and 
250 km s“^. To investigate our sensitivity to different 
parts of the sky, we cut our sample in half in longitude 
and in latitude. 

We also investigate different model specifications, and 
compare the results from the restricted and unrestricted 
model for each sample. Finally, the location and size of 
the bins used in the estimations may affect the results. 
We therefore re-estimate the model parameters using an 
Epanechnikov kernel centered at R with finite support 
on a total bandwidth of 4 kpc chosen to match the size 
of fixed bins used. This method gives greater weight to 
stars near R, rather than weighting them uniformly as 
occurs with fixed bins, and provides continuous estimates 
of the parameters as functions of R. 

In every case, the resulting parameter estimates change 
by less than about one standard deviation from our orig¬ 
inal results with the exception of a systematic distance 
error, where the change in parameter estimates is less 
than two standard deviations. We conclude that our re¬ 
sults are numerically robust within the estimated errors. 


TABLE 1 
Tilt Angles 


R 

N 

“He 

<^R4> 

aecjy 

(kpc) 


(deg) 

(deg) 

(deg) 

16.7 

2482 

Q q+2.9 

Z.O_2.9 

1 q+5.9 
-*-•^-5.9 

6.51S;? 

22.4 

2483 

—n 1+^-^ 

—0 

36.3llt:6 


gles, Oy, derived bv iSmith et al.l (1200911 . 

2E- ■ 

tan {2aij) = , (4) 

^jj 


where is the covariance between the velocity compo¬ 
nents Vi and Vj. The tilt angles specify the orientation 
of the velocity ellipsoid as the angle between the i-axis 
and the major axis of the ellipse resulting from the pro¬ 
jection of the three di mensional velocity ell i psoid onto 
the ij -plan e (see, e. o.. iBinnev &: Merrifieldl (119981 1 and 


iSmitir et al.l (|2009ll l. iSmith et al.l ( 200911 showed that if 
the inner halo is in steady state and the velocity ellipsoid 
is everywhere aligned with the Galactic spherical coor¬ 
dinate system, then the gravitational potential must be 
spherically symmetric. This result holds for the velocity 
ellipsoid of any tracer population, whether its density 
distribution is oblate, prolate or triaxial. 

We find no evidence of any clear tilt in our analysis 
of the combined sample using the eight disjoint, approx¬ 
imately equally populated bins of 2,842 or 2,843 stars 
in the interval 6 < i? < 30 kpc. None of the esti¬ 
mated tilt angles differs from zero at a Icr statistical 
significance level. For R < 15 kpc, the tilt angles are 
poorly estimated: the Icr confidence intervals for the es¬ 
timates of the tilt angles are of order of tens of degrees. 
Our most tightly constrained estimates occur for_the bins 
with mean R = 16.7 (R G (15.3,18.5)) kpc and R = 22.4 
{R G (18.5, 30.0)) kpc. The corresponding estimates for 
the tilt angles with Icr errors are given in Table [TJ None 
is significantly different from zero even at Icr. The tilt 
angle aetf, is less well constrained due primarily to the 
large uncertainty in estimated covariance between the 
two components. 

Our findings are consistent with a spherically sym¬ 
metric gravitational potential in t he inner ha lo, _ as 

suggested by several recent studies (ISmith et al.l 120091 : 


iKoposov et al.l[2?)Tol : lAgnello fc Evansl2012ll . Ourresults 

accord with those of ISmith et al.l (|2009fl . who found tilt 
angles consistent with spherical symmetry for ^ 1,500 
nearby halo subdwarf stars with heliocentric distances of 
< 5 kpc (and A < 11 kpc) along the ^ 250 deg^ cov¬ 
ered by SDSS Strip 82. Their work was limited both 
in distance and sky coverage. Our sample, on the other 
hand, extends to A = 30 kpc with substantially greater 
sky coverage, including high Galactic latitudes where any 
contribution from the Galactic disk is negligible. 


4.5. Alignment of the Halo Velocity Ellipsoid and the 
Gravitational Potential of the Milky Way 

The alignment of the velocity ellipsoid for halo stars 
with respect to the Galactic coordinate system provides 
a powerful probe of the gravitational potential of the 
Galaxy. The alignment can be described by the tilt an- 


4.6. Measuring the Interior Mass of the Milky Way 

Many studies have exploited the dyna mics of halo stars 
to measure the mass of the Milky Wa y (jXue et al.ll2008l : 
iDeason et al.l 120121 : iKafle et all 1201211 by assuming dy¬ 
namical equilibrium and fitting specific forms of the dis¬ 
tribution function, postulated to depend only on two in¬ 
tegrals of motion, the binding energy and the angular 
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momentum. The mass distribution of the Milky Way is 
commonly measured using the steady-state Jeans equa¬ 
tion for a spherical potential usi ng a tracer population of 
stars (|Binnev fc Tremain3l2n08l eq. 4.37): 


M{< R) 


( d In ptr 
dlni? 


din (T^ 
dlni? 



(5) 


The mass interior to a radius i? is a function of the 
anisotropy, j3, and the logarithmic gradients of the ra¬ 
dial velocity dispersion, an, and the density of trac- 


has been extensively studied {e.q-< 

Yannv et al.l 120001: 

Chen et al.l 120011: INewbersf & Yannv 

120061: lJuric et al.l 

20081: iDeason et al.ll20111 leading to a recent consensus 


120Ilf) for the inner halo, R <27 kpc. 


p{Rq)(xRq°', ( 6 ) 

i? 2 = + + 

r 


where the power law index is a = —2.3 to —2.6 and the 
minor axis to major axis ratio is g = 0.6 — 0.7. Recent 
simulations by IWang et al.l (I2015f) evaluate potential bi¬ 
ases in estimating the mass of the Milky Way using dy¬ 
namical tracers. They find that although deviations from 
spherical symmetry are relatively unimportant, devia¬ 
tions from dynamical equilibrium can cause significant 
bias. To avoid any discontinuity in the break region, 
where the assumption of dynamic equilibrium may not 
hold, we calculate the mass within 12 kpc of the Galactic 
center, M {R < 12 kpc), using the results from the 4 kpc 
wide bin centered on 12 kpc and the Jeans equation ([5]). 
We find M(i? < 12 kpc) = 1.3 x 10“ Mq. 

We estimate that the mass uncertainty is 40%, arising 
largely from the uncertainty in the dlncr^/dlni? term in 
Equation [51 The estimated error in is only 10%, and 
reported values of d In ptr/d In R also vary by 10%. When 
we compare adjacent bins to i? = 12 kpc to estimate the 
scatter in dlncr^/dlni?, however, we find that this term 
varies by 35%. To compare with previous work, we also 
calculate the mass for the R = 25 kpc bin. We find 
M {R<25 kpc) 2.6 x 10^^ Mf?^, a result that agrees to 
within 20% with iKafie et al.l (l2012fl . 


5. DISCUSSION 

Motivated by the large tangential velocity dispersion 
estimates in the region 15 ;< i? < 25 kpc, we exam¬ 
ine whether the observational data are consistent with 
our one fundamental assumption: that the three compo¬ 
nents of velocity, (u/{,ue,u$), are normally distributed. 
The normal distribution is described by two parameters, 
its mean and dispersion. If the distribution of halo star 
velocities is not intrinsically normal and we (incorrectly) 
try to represent it as such, we expect the mean will re¬ 
main largely unchanged but the dispersion will vary. No¬ 
tably, the Sgr stream passes through our survey region 
in a way that introduces a bimodal velocity distribution 
of stars into the sample beginning at i? ^ 15 kpc. 


5.1. Departures from Normality 

If the underlying theoretical assumptions are correct, 
the observed line-of-sight velocities should be normally 



Vrf [km s-'] 

Fig. 13.— Top panel: Observed Vj.f distribution of stars in 
the bin 14 < i? < 18 kpc, superi mposed with a fitted , three- 
component, normal mixture model IIBenaglia et al.l 1200911 . There 
appears to be an excess of stars in both the negative and posi¬ 
tive v elocity wings of the distr ibution. Bottom panel: Distribution 
of the ILaw fc Maiewskil <1201111 1 Sgr A^-body model sampled in the 
identical way for 14 < ij < 18 kpc. The simulated Sgr stars have 
a bi-modal distribution of line-of-sight velocities in this region. 

distributed since linear combinations of normally dis¬ 
tributed random variables are also normally distributed. 
Yet Section 12.41 shows that the distribution of observed 
radial velocities departs from normal in certain regions 
and fails standard statistical tests for normality. Physi¬ 
cally, known or unknown star streams or other velocity 
structure may be responsible for this observed departure 
from normality. Correlations between angular position 
and velocity break our model assumptions. 

To better understand the possible causes of these dis¬ 
crepancies, we look more closely at the break region by 
selecting the 4,077 stars in the 14 < i? < 18 kpc bin 
and plotting the histogram of Vrf in Figure [T31 The 
histogram reveals deviations from a normal distribution 
both in the center and in the wings of the distribution. 
The Kolmogorov-Smirnov test rejects normality at a rea¬ 
sonable level of statistical significance (p = 0.069), and 
the more powerful Anderson-Darling test that gives more 
weight to the tails than the Kolmogorov-Smirnoff test 
strongly rejects normality (p = 9.9 x 10“®). 

5.2. Non-Normality of Sgr Stream Velocities 

At least one significant halo structure lies within our 
survey footprint: the Sgr stream. To assess the im- 
pa ct of the Sgr stream, w e turn to the iV-body model 
of ILaw fc Maiewskil ([201011 . Their model is designed to 
match existing observational constraints on the location 
and motion of the Sgr stream in a triaxial potential. We 
sample the fV-body model for test particles that he at 
\Z\ > 5 kpc in the range 14 < i? < 18 kpc and that 
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fall within an approximation of our survey footprint on 
the sky (6 > 30° or (6 < —30° and I G (50°, 200°))). 
Although the majority of Sgr is at larger distances, a 
fraction of the stream is present in the 14 < ii < 18 kpc 
region. Figure [T3] (lower panel) shows the Vrf distribu¬ 
tion of Sgr iV-body particles in this region. 

Because the Sgr stream wraps around the sky in a 
roughly polar orbit, the Sgr stars passing through the 
14 < A < 18 kpc region exhibit multiple velocity peaks 
and a broad velocity dispersion. Suggestively, the two 
main peaks in Vrf from the Sgr model coincide with the 
excesses found in the wings of the observed halo star dis¬ 
tribution (Figure [T5) in a fitted, three co mponent, nor¬ 
mal mixture model ( Benaglia et al.l[2009ll . 

Previous research estimates that up to 10% of the stars 
in the survey reg ion may originate in the Sgr stream 
((King et al.ll2012ri . The presence of Sgr stream stars in 
our sample may thus account, at least partially, for the 
observed discrepancies in the observed velocity disper¬ 
sions and anisotropy of the halo in the break region. We 
test this hypothesis by removing all stars in our sample 
within 10° of the Sgr stream, \B\ < 10 °, using the Sgr 
stream coordinates, {A,B), defined bv iBelokurov et ahl 

(IMl . 

Removing the region of the Sgr stream from our sample 
modestly lowers our estimates of the dispersions, which, 
in turn, increases our estimate of the anisotropy. The 
anisotropy in the 14 < i? < 18 kpc bin changes from 
/3 = —8 to P = —6, in better agreement with previ¬ 
ous results (Figure IT^ . The change in P is formally 
less than Itr, but there may be other star s tream s for 
which we do not accoun t. Indeed, iBell et aH (|2008ll and 
iSchlaufman et al.1 (1200911 identify 30% - 40% of F-type 
stars at 15 kp c depths in c o heren t spatial or velocity 
structures, and l.lanesh et al.1 (|2015l l show that the per¬ 
centage of halo stars in coherent structures increases with 
depth. We conclude that the traditionally assumed nor¬ 
mal velocity distribution model may not properly repre¬ 
sent the substructure of the stellar halo. 

6. CONCLUSION 

We use a large sample of 19,859 stars at 6 < i? < 30 
kpc to investigate the mean velocities, velocity disper¬ 
sions, covariances, and anisotropy of the Milky Way halo. 
This dense sample enables finer binning than previously 
used and allows us to investigate the 15 < i? < 25 kpc 
anomaly in anisotropy at higher statistical significance. 

We begin by presenting a new radial velocity survey 
of 6,174 faint F-type stars observed with the Hectospec 
spectrograph using the MMT telescope. F-type stars are 
dense tracers of both the thick disk and halo. To focus 
on halo kinematics, we restrict our analysis to stars with 
\Z\ > 5 kpc. We add stars from published SDSS radial 
velocity samples to create a combined sample of 19,859 
stars that span 6 < i? < 30 kpc. 

We use the Sun’s offset from the Galactic center to re¬ 
cover tangential velocity information from the observed 
line-of-sight velocities utilizing standard statistical meth¬ 
ods. We make the minimal assumption that the under¬ 
lying stellar velocity distribution is normal. We use a 
maximum likelihood procedure to calculate the veloc¬ 
ity means, dispersions, covariances, and anisotropy. We 
find that the alignment of the velocity ellipsoid is consis¬ 
tent with a spherically symmetric gravitational potential. 


From the spherical Jeans equation, we estimate the mass 
of the Milky Way within 12 kpc is M {R < 12 kpc) = 
1.3 X 10^^ Mq with an uncertainty of 40%. 

A significant region of discontinuity 15 < i? ;< 25 
kpc exists where the estimated velocity dispersions and 
anisotropy diverge from their anticipated values, confirm¬ 
ing the break region observed by others. The estimated 
tangential velocity dispersions in this region are so large 
that stars would be unbound, an unphysical result. Yet 
the results are numerically robust. In sensitivity analy¬ 
ses (ie., using a different solar motion, different distance 
scale, different survey footprint, etc.), the maximum like¬ 
lihood calculation yields parameters that change by less 
than about one standard deviation from our original re¬ 
sult. Conversely, if we input simulated data drawn from 
known normal velocity distributions, the maximum likeli¬ 
hood estimation finds the correct velocity dispersion and 
anisotropy parameters. 

We suggest that the discontinuity in the region 15 < 
R < 25 kpc arises from the failure of the normal distri¬ 
bution model to describe the actual velocity data. Phys¬ 
ically, known or unknown star streams or other velocity 
substructure may be responsible for the departure from 
normality. The predicted contribution of the (polar or¬ 
biting) Sgr stream in our survey region, for example, is 
a bi-modal distribution of stars in the wings of the ob¬ 
served radial velocity distribution. Sgr by itself cannot 
explain the discontinuity, but Sgr is unlikely to be the 
only structure in the halo. 

The upshot is that larger radial velocity samples alone 
cannot improve our understanding of the halo using the 
standard statistical approach. Significant improvement 
requires direct tangential velocity constraints for halo 
stars, like those soon to be provided by Gaia. The com¬ 
bination of our radial velocity measurements with Gaia 
proper motions will thus be very useful for understanding 
the physical nature of the 15 < i? < 25 kpc discontinuity 
region and the kinematics of the Milky Way halo. 
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APPENDIX 

APPENDIX 

Coordinate Systems and Projection Factors 

We employ four coordinates systems: heliocentric Cartesian {x, y, z) and spherical (d, b, 1) coordinates, and Galacto- 
centric Cartesian {X, Y,X) and spherical {R, 6, (j) coordinates. In the Cartesian coordinate systems, the Galactocentric 
X- and V-axes lie in the Galactic plane; the heliocentric xy-plane is parallel to, but slightly above, the Galactic plane 
since the Sun lies slightly above it. The z- and Z-axes are perpendicular to the xy- and W-planes, respectively, 
forming right-handed coordinate systems. The positive x- and V-directions are defined as from the Sun toward the 
Galactic cen ter. Since the Sun is located 8 kpc from the Galactic center and 0.0196 kpc above the Galactic plane 
(|Reedl[200^ . the two Cartesian coordinate systems transform as 

X = x — 8.0, Y = y, Z = z-|-0.0196 kpc (Al) 

The spherical coordinates are defined in terms of their Cartesian coordinates with radii d = ^x'^ + y"^ + z"^ and 
R = y/X"^ -I- _|_ ^2 fQj. |-]^g heliocentric and Galactocentric systems, respectively. Longitudes I and (j) are measured 

in the xy- or XY- plane from the x- or X-axis toward the y- or V-axis. The second angle of the spherical coordinates 

is defined differently for the heliocentric and Galactocentric systems. For the heliocentric system, b is dehned as 
the latitude, the angle measured from the ly-plane with the direction toward the z-axis taken as positive. For the 
Galactocentric system, 0 is dehned as the colatitude, the angle measured from the Z-axis. 

The line-of-sight velocity component projection factors are: 

Ph = cos b cos I sin 9 cos 4> -\- cos b sin I sin 9 sin (j) -|- sin b cos 9 (A2) 

pe = cos b cos I cos 9 cos (f) -\- cos b sin I cos 9 sin (j — sin b sin 9 (A3) 

p^ = — cos b cos lsm(j)-\- cos b sin I cos (j) (A4) 

A projection factor pg = 0.3, for example, means that a tangential velocity component of vg = 100 km s“^ contributes 
30 km s“^ to the line-of-sight velocity, vios- 


Probability Density Function 

Following standard theoretical development, we assume that the components of stellar velocity, Vi, are normally 
distributed with means, pi, and standard deviations, Ui, 


ViN , iG{R,9,(j)}. (A5) 

The observed line-of-sight velocity for a given star is then the sum of the projections of each velocity component 
onto the line of sight, as previously dehned in Equation [5] Because the line-of-sight velocity is a linear combination 
of the velocity components, assumed to be normally distributed, the observed line-of-sight velocity is also normally 
distributed with mean, 

E[vios] = PrPr + pope + = pV. (A6) 


and variance, 

Var[vios\ = Ear[p'v] = p'Sp, (A7) 

where is the variance-covariance matrix among the Vi and is symmetric, = Yji. Gombining these results 

yields the distribution of line-of-sight velocities 

vios ~ N (p'p, p'Sp) (A8) 

= N{prPr -I- pepig -I- P4,p^, 

P^rYrr -\- p^gYgg -\- -|- 2pRpgT,Rg -\- 2 prP^Yr^ -\- 2pgp^T,g^) 


Now consider the probability / (a:|0) for hnding a star at Galactic coordinates (i?, 9, (j>) with line-of-sight velocity vios- 
Dehning an observation as x = {R,0,ji,vios} and the model parameters as 0 = {p,S}, the probability distribution 
function is 


f(xl&) 



v^Vp'Sp 


exp 


1 (vios - P'/J-f 

1 

1 (vios - p'pf 

2 p'Sp 


2 p'Sp 


(A9) 


where p ^R^ is the number density of stars. This generalization of the Schwarzschild distribution, originally proposed 
bv ISchwarzschildl (|1907f) to model the velocity distribution of stars in the solar neighborhood, further allows correlations 
among the spatial velocities. Since p ^rJ is not a function of the model parameters, it may be omitted as it does 
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TABLE 2 

Hectospec F Stars 


R.A. 

(h:m:s) 

Decl. 

(d:m:s) 

^helio 

(km s“^) 

ro 

(mag) 

Mr 

(mag) 

R 

(kpc) 

Z 

(kpc) 

0:39:30.452 

-1-25:14:52.60 

-287.6 ± 19.8 

19.940 ± 0.027 

5.241 

13.94 

-5.31 

0:39:59.518 

-1-25:12:41.95 

-44.8 ± 18.7 

20.182 ±0.028 

5.821 

12.89 

-4.54 

0:41:30.298 

-1-25:09:09.98 

-167.8 ±26.6 

20.520 ± 0.043 

6.412 

12.25 

-4.05 

0:41:43.272 

-1-25:55:22.42 

-224.6 ± 12.4 

19.196 ±0.018 

4.790 

13.07 

-4.57 

0:41:44.054 

-1-25:01:11.58 

-75.9 ±23.3 

20.510 ± 0.037 

5.247 

16.22 

-6.92 


Note. — (This table is available in its entirety in machine-readable and Virtual 
Observatory forms in the online journal. A portion is shown here for guidance regarding 
its form and content.) 


not enter into the subsequent maximum likelihood analysis, an advantage of the maximum likelihood technique. We 
estimate nine parameters of 0 in our statistical model: three means /ii and six independent elements of the symmetric 
covariance matrix S. 

Denoting an individual observation by a; = {R,9,(p,vios} and a set of observations by x = {xi,X 2 ---Xn}, the log 
likelihood function is 

n 

C{&\x) = '^logf {x^\&) . (AlO) 

Data Table 

Table m presents the clean sample of 6,174 F-type stars from the Hectospec radial velocity survey. For each star 
we list the position (epoch J2000), our heliocentric radial v elocity measuremen t, the SDSS DRIO de-reddened r-band 
magnitude, our estimated absolute magnitude Mr using the llvezic et al.l (I2008D photometric parallax relation, and the 
corresponding Galactocentric radial R and vertical Z distances, calculated assuming the Sun is at i? = 8 kpc. The full 
version of Table [5] is available in the online journal. 


Results Table 

Table [3] presents the estimation results for the Hectospec F star and combined samples. 
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TABLE 3 

Estimation Results 


N 

Rmin 

(kpc) 

Rmax 

(kpc) 

R 

(kpc) 

(km s“^) 

(km s~^) 

(^(f) 

(km s“^) 

(km s“^) 

(km s“^) 

^ <p 

(km s“^) 

(km^ s“^) 

(km^ s“^) 

(km2 s-2) 

/3 


Hectospec F Star Sample 

762 

6.0 

10.6 

9 .3 




130.3 

111.5 

90.0 




0.40 








(20.6) 

(10.7) 

(29.5) 




(0.38) 

762 

10.6 

12.4 

11.4 




157.5 

42.5 

72.8 




0.86 








(21.6) 

(93.2) 

(58.2) 




(0.33) 

762 

12.4 

14.9 

13.5 




94.7 

167.0 

165.6 




-2.08 








(26.7) 

(41.4) 

(29.6) 




(2.91) 

763 

14.9 

29.8 

18.2 




98.9 

203.3 

132.3 




-2.01 








(9.2) 

(30.7) 

(35.0) 




(1.44) 

Combined Sample - equally populated bins 

2483 

6.0 

9.6 

8.4 

0.8 

-4.9 

-7.0 

155.3 

109.8 

88.3 

428.7 

-271.8 

-80.4 

0.59 





(9.1) 

(5.6) 

(7.1) 

(14.0) 

(7.8) 

(18.6) 

(1971.1) 

(2251.1) 

(1670.5) 

(0.13) 

2482 

9.6 

10.6 

10.1 

-4.6 

-8.0 

3.9 

156.5 

98.2 

86.2 

-123.2 

2442.2 

-456.0 

0.65 





(6.9) 

(5.7) 

(7.4) 

(16.5) 

(14.2) 

(28.7) 

(1949.6) 

(2944.3) 

(2712.9) 

(0.19) 

2482 

10.6 

11.5 

11.1 

-8.5 

-4.7 

-7.2 

107.4 

117.6 

110.3 

-2806.8 

-2923.5 

-6839.6 

-0.13 





(5.1) 

(5.4) 

(7.7) 

(21.5) 

(20.3) 

(26.3) 

(1190.2) 

(2357.3) 

(2723.3) 

(0.83) 

2482 

11.5 

12.5 

12.0 

4.9 

3.5 

3.8 

150.1 

38.5 

85.9 

-2088.1 

530.2 

-4335.6 

0.80 





(4.6) 

(5.7) 

(9.0) 

(12.3) 

(76.6) 

(40.2) 

(1071.4) 

(2054.6) 

(2761.1) 

(0.29) 

2483 

12.5 

13.7 

13.1 

-8.1 

-1.2 

-8.6 

105.0 

165.4 

194.0 

1896.4 

5801.2 

5498.6 

-1.95 





(3.9) 

(6.1) 

(9.7) 

(11.8) 

(17.8) 

(17.9) 

(913.9) 

(1961.8) 

(3269.4) 

(1.11) 

2482 

13.7 

15.3 

14.4 

-8.3 

6.3 

-12.7 

95.1 

205.3 

172.6 

3053.6 

1448.5 

2843.3 

-2.97 





(3.6) 

(7.3) 

(9.4) 

(10.6) 

(18.3) 

(19.2) 

(1023.6) 

(1732.2) 

(4172.5) 

(1.52) 

2482 

15.3 

18.5 

16.7 

-8.6 

-19.4 

-40.1 

56.8 

256.0 

225.8 

1243.5 

494.8 

831.7 

-17.08 





(2.9) 

(8.3) 

(9.5) 

(9.3) 

(14.4) 

(14.2) 

(793.6) 

(1257.0) 

(4391.9) 

(7.50) 

2483 

18.5 

30.0 

22.4 

-10.7 

-24.8 

-56.6 

76.8 

259.1 

195.8 

57.2 

30.9 

10559.0 

-7.95 





(2.4) 

(11.3) 

(10.1) 

(4.7) 

(18.8) 

(21.7) 

(952.3) 

(1046.3) 

(6406.6) 

(2.32) 







Combined Sample 

- 4 kpc bins 






3356 

6 

10 

8.8 

-7.4 

-8.6 

-9.9 

149.4 

109.4 

91.5 

65.0 

-438.2 

-212.8 

0.54 





(7.3) 

(4.8) 

(6.1) 

(11.6) 

(6.6) 

(15.8) 

(1566.3) 

(1951.8) 

(1512.3) 

(0.12) 

5680 

7 

11 

9.7 

-3.6 

-6.3 

-4.7 

137.7 

109.8 

110.5 

-505.9 

1309.8 

-1025.6 

0.36 





(4.6) 

(3.6) 

(4.7) 

(7.5) 

(5.3) 

(11.0) 

(987.4) 

(1460.6) 

(1263.4) 

(0.12) 

7915 

8 

12 

10.4 

-2.3 

-4.1 

-2.9 

129.6 

109.9 

112.6 

-1373.8 

447.5 

-2247.6 

0.26 





(3.3) 

(3.0) 

(4.2) 

(5.2) 

(4.8) 

(9.4) 

(692.5) 

(1124.0) 

(1088.7) 

(0.11) 

9505 

9 

13 

11.1 

-2.8 

-1.8 

-4.2 

128.6 

114.8 

123.9 

-603.9 

1515.5 

-1374.5 

0.14 





(2.7) 

(2.8) 

(4.0) 

(3.9) 

(4.9) 

(8.6) 

(594.3) 

(1023.9) 

(1112.6) 

(0.11) 

9658 

10 

14 

11.9 

-3.8 

-2.4 

-2.1 

129.3 

112.7 

125.0 

-426.5 

1837.7 

-1941.9 

0.15 





(2.3) 

(2.8) 

(4.2) 

(3.5) 

(6.2) 

(9.5) 

(518.2) 

(955.3) 

(1186.2) 

(0.12) 

8606 

11 

15 

12.8 

-6.2 

-2.2 

-6.3 

123.5 

126.3 

135.3 

315.9 

1725.8 

-471.1 

-0.12 





(2.2) 

(3.2) 

(4.8) 

(3.9) 

(8.2) 

(11.0) 

(533.1) 

(989.9) 

(1543.6) 

(0.19) 

6933 

12 

16 

13.7 

-6.1 

3.8 

-11.8 

110.4 

157.7 

166.7 

1687.6 

2648.4 

1922.0 

-1.16 





(2.2) 

(3.9) 

(5.6) 

(4.8) 

(9.6) 

(11.0) 

(578.4) 

(1082.8) 

(2092.4) 

(0.39) 

5354 

13 

17 

14.6 

-6.2 

-0.7 

-12.3 

97.2 

190.1 

180.4 

2195.8 

2557.4 

3649.3 

-2.64 





(2.3) 

(4.8) 

(6.5) 

(5.8) 

(11.1) 

(12.4) 

(641.7) 

(1142.6) 

(2718.3) 

(0.79) 

4077 

14 

18 

15.6 

-5.9 

-3.1 

-21.1 

72.0 

233.9 

201.0 

1907.1 

1291.7 

2263.5 

-8.16 





(2.4) 

(6.0) 

(7.4) 

(6.9) 

(10.7) 

(12.3) 

(696.4) 

(1125.7) 

(3316.7) 

(2.46) 

3069 

15 

19 

16.7 

-8.3 

-10.4 

-39.6 

56.5 

257.9 

215.8 

1672.3 

-556.6 

-1410.4 

-16.71 





(2.6) 

(7.5) 

(8.4) 

(8.3) 

(12.5) 

(12.9) 

(711.3) 

(1107.7) 

(3941.0) 

(6.56) 

2504 

16 

20 

17.7 

-13.3 

-29.2 

-50.8 

53.1 

270.2 

217.6 

1050.0 

-1129.9 

-84.7 

-20.35 





(2.6) 

(8.9) 

(9.2) 

(7.2) 

(12.8) 

(13.4) 

(750.1) 

(1058.7) 

(4577.8) 

(7.40) 

2013 

17 

21 

18.7 

-12.7 

-25.2 

-59.2 

58.5 

257.8 

218.5 

334.5 

-1897.6 

-1140.2 

-15.68 





(2.8) 

(10.2) 

(10.2) 

(6.8) 

(16.5) 

(15.7) 

(786.5) 

(1093.2) 

(5339.4) 

(5.60) 

1615 

18 

22 

19.7 

-15.7 

-23.7 

-63.7 

58.9 

284.7 

217.0 

235.7 

-1582.4 

4292.7 

-17.48 





(3.1) 

(12.7) 

(11.7) 

(7.4) 

(19.3) 

(19.4) 

(986.6) 

(1276.8) 

(7048.2) 

(6.75) 

1328 

19 

23 

20.7 

-14.3 

-37.3 

-65.4 

52.6 

298.9 

243.7 

-1262.7 

-125.2 

10052.2 

-25.89 





(3.2) 

(14.3) 

(13.4) 

(8.5) 

(21.4) 

(21.6) 

(1063.2) 

(1321.9) 

(7893.0) 

(11.88) 

1093 

20 

24 

21.8 

-8.9 

-12.4 

-62.1 

54.4 

323.9 

241.3 

66.4 

-920.7 

3538.6 

-26.59 





(3.6) 

(17.1) 

(15.2) 

(9.9) 

(26.8) 

(25.2) 

(1355.4) 

(1431.5) 

(8801.9) 

(13.94) 

918 

21 

25 

22.8 

-9.9 

-26.5 

-75.5 

48.7 

362.1 

248.9 

1689.7 

462.6 

13621.3 

-39.73 





(4.0) 

(20.2) 

(16.7) 

(11.8) 

(28.1) 

(26.8) 

(1763.6) 

(1537.5) 

(10629.8) 

(25.11) 

763 

22 

26 

23.7 

-2.2 

-21.5 

-48.2 

62.8 

329.5 

232.4 

959.3 

550.4 

7405.9 

-19.61 





(4.3) 

(22.2) 

(18.4) 

(10.4) 

(31.7) 

(34.1) 

(1954.0) 

(1713.7) 

(12355.0) 

(10.39) 

649 

23 

27 

24.8 

-0.9 

-12.6 

-50.3 

85.3 

279.1 

142.4 

1262.5 

873.7 

19147.5 

-5.75 





(4.8) 

(25.4) 

(20.7) 

(10.1) 

(46.1) 

(70.0) 

(2318.5) 

(1984.0) 

(15005.4) 

(4.26) 

538 

24 

28 

25.8 

-1.5 

-46.9 

-28.4 

92.6 

191.3 

123.6 

-1807.9 

2175.8 

26741.0 

-2.02 





(5.0) 

(26.9) 

(23.0) 

(11.4) 

(74.5) 

(112.0) 

(2260.6) 

(2252.1) 

(16965.0) 

(3.70) 

443 

25 

29 

26.7 

-6.1 

-36.1 

-11.1 

90.4 

140.5 

160.5 

-5339.9 

-251.5 

-2283.8 

-1.78 





(5.5) 

(30.3) 

(28.0) 

(14.1) 

(138.9) 

(114.7) 

(2179.6) 

(2674.5) 

(19133.2) 

(5.12) 












